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We present high precision results on the static qq potential on 32 4 and smaller lattices, using the standard 
Wilson action at /3 = 6.0, 6.2, 6.4, and 6.8 on the Connection Machine CM-2. Within our statistical errors (~ 1%) 
we did not observe any finite size effects affecting the potential values, on varying the spatial lattice extent from 
0.9 fm up to 3.3 fm. We find violations of asymptotic scaling in the bare coupling up to f3 = 6.8. We demonstrate 
that scaling violations on the string tension can be considerably reduced by introducing effective coupling schemes, 
which allow for a safer extrapolation of A_l to its continuum value. We are also able to see and to quantify the 
running of the coupling from the interquark force. From this we extract the ratio y^/Az, ■ Both methods yield 



consistent values for the A-parameter: A-rrc = 0.558 



+0.017 



246+lMeV. 
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1. Introduction 

It is of great interest to extend lattice tech- 
niques into a regime where perturbative meth- 
ods can be applied |I|-^J]. We achieve a lattice 
spacing as small as a m (6 GeV)" 1 on a 32 4 lat- 
tice at p = 6.8. From the static ST/ (3) quark- 
antiquark potential we can extract the running 
coupling "constant" within an accuracy that com- 
petes with current (real QCD) experiments. 

2. Simulation 



In our simulations [|5j,|6| we have realized lattice 
volumes of L| x T T = 32 4 for the /?- values 6.0, 
6.4, and 6.8. and a 24 3 x 32 lattice for /3 = 6.2. 
The corresponding spatial lattice extents in order 
of increasing /3 are: 3.25, 2.74, 1.74, and 1.05 fm 
with the physical scale set by the value y/a = 440 
MeV for the string tension. In order to investigate 
finite size effects (FSE) we simulated, moreover, 
a 16 4 lattice at P = 6.0, 16 3 x 32, 24 3 x 32, and 
32 3 x 16 lattices at p = 6.4, and a 16 3 x 64 lattice 
at P = 6.8. 

We used the hybrid overrelaxation algorithm 
for updating the gauge fields, as advocated in 
Kennedy's review talk M. Great care was taken 
in optimizing the overrelaxation step || . 

To improve the projection of our operators onto 



the qq ground state we applied an iterative lo- 
cal gauge covariant smoothing procedure || on 
the spatial links before constructing Wilson loops. 
Overlaps of 95(80)% for small (large) spatial qq 
separations R were reached. Various smoothened 
on- and off-axis Wilson loops were measured ev- 
ery 100 sweeps. This separation was found to 
be sufficient to guarantee an (almost) indepen- 
dent sampling. Remaining autocorrelation effects 
were healed by blocking the data into bins of rea- 
sonable length before statistical analysis. On all 
volumes 0(100) measurements were performed. 

The spatial qq separations R = ne, with ej = 
(1, 0, 0), (1, 1, 0), (2, 1, 0), (1, 1, 1), (2, 1, 1), (2, 2, 1) 
were realized, n was increased up to Ts/2 for 
i = 1,2,4, and up to Lg/4 for the remaining di- 
rections. Altogether this yields 72 different values 
of R on the 32 3 x Lt lattices. The time separa- 
tions T — 1, 2, .... 10 were used to test exponen- 
tial behaviour and for extraction of the potential 
values. 

3. Results 

qq potential: We follow Chris Michael [l]] and 
start from the ansatz 

V(R) = V Q +KR-e (±j+ + IG L (R)\ +-^.(1) 
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The lattice propagator for the one gluon exchange 




Figure 1. Potential V C (R) at (3 = 6.4, V = 32 4 



Table 1 

Fit parameters (see Eq. (1)). 





= 6.0 


/3 = 6.4 


K 


0.0513(25) 


0.0148(3) 


e 


0.275(28) 


0.315(15) 


V 


0.636(10) 


0.601(4) 


I 


0.64(12) 


0.56(6) 


f 


0.041(58) 


0.075(18) 




/3 = 6.5[2J 


13 = 6.8 


K 


0.0114(2) 


0.0053(2) 


e 


0.311(14) 


0.311(10) 


V 


— 


0.549(2) 


I 


0.64(6) 


0.56(4) 


f 


0.067(13) 


0.094(13) 



first two data points are excluded. 

We have included the corresponding fit param- 
eters from the recent UKQCD ||] investigation on 
a 36 4 lattice at /3 = 6.5 in Table 0. For (3 > 6.4 
all Coulomb coefficients e are definitely different 
(but remarkably stable) from the string vibration 
value 7r/12 « 0.262. Vo decreases with the cou- 
pling. We emphasize that the parameter / is es- 
tablished to increase with (3 as expected. I seems 
to vary only slightly. 

Finite size effects: On varying the lattice ex- 
tent at constant lattice spacing from 1.7 to 3.3 fm 
(/3 = 6.0), and from 0.9 up to 1.8 fm (/3 = 6.4) 
we found no evidence of FSE on our potential 
data within statistical accuracy (« 1%). How- 
ever, a comparison of results on the (3 = 6.8 lat- 
tices (0.53 and 1.05 fm) shows deviations (up to 
3cr), mainly for small R separations where statis- 
tical errors are small. We conclude that a lattice 
extent aLs > 0.9 fm fw 2/ v / cr suffices in our case. 
The same limit was also found in Ref. ||. This 
condition coincides with the transition into the 
deconfined phase at clLt > 1/Tc ~ 2/ ' yja. 

Note, that a small volume can also cause prob- 
lems of another kind: It is difficult to fit the 
potential to a small number of data points with 
some confidence. This effect becomes dominant 
for (3 < 5.9 where the physical limit is L$ > 8. 
In any case, it pays to work on a large lattice. 
The computational task (except memory require- 
ment) remains comparable because self averaging 
suppresses fluctuations and reduces the number 
of required Monte Carlo sweeps. 



Gl(R) has been calculated in the large volume 
limit. The parameter I controls violations of ro- 
tational symmetry within this ansatz. The term 
f/R simulates deviations from a pure Coulomb 
behaviour and is expected to be positive from 
asymptotic freedom, and to increase with the en- 
ergy resolution. 

In Figure 1 the lattice corrected data points 
Vc(R) = V(R) + el(G L (R) - l/R) are plotted 
together with the interpolating fit curve. The fig- 
ure demonstrates the success of ansatz Eq. (1) as 
all data points line up on a smooth curve. Our 
potential fits yield x 2 /Nof < 1 as long as the 



Table 2 

in units of the string tension. 



AZ 1 



(3 Ref y/B/h L v^MT 



5.7[12] 

5.8[12] 

5.9g| 

6.0 

6.2 

6.4 

6.5 || 

6.8 



124.7 (07) 
112.4 (10) 
102.9 (14) 
96.5 (23) 
86.4 (10) 
81.3 (08) 
80.0 (14) 
76.9 (13) 
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63.3 (04) 
63.0 (06) 
61.2 (08) 
60.0 (15) 
56.9 (07) 
55.7 (05) 
55.7 (10) 
55.7 (09) 



53.2 



7.3 



^JW 1 



55.7 (03) 
55.6 (05) 

54.3 (07) 

53.4 (13) 
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50.1 (09) 
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Figure 2. Continuum extrapolations of A L . 



Scaling: Within our range of lattice spacings, 
scaling of physical quantities works nicely: In the 
present investigation all potential data, measured 
at different (3 values, fall onto a single curve when 
scaled to Vb = 0, and K = 1. However, though 
scaling seems to be restored, it still does not 
follow the perturbatively expected two-loop be- 
haviour a a = /(/3)/Ai with a a = \f~K j 'y/a, and 
f((3) — the integrated two-loop (5- function. This 
is reflected in the first column of Table 0. It was 
argued that the lack of asymptotic acaling may be 
caused by a bad choice of expansion parameter. 
For remedy, a variety of schemes with redefined 
expansion parameters have been proposed in the 
past |,0. 

Here, we attempt to achieve an improvement 
by using the (inverse) coupling (3e that has origi- 
nally been invented by Parisi |9| , and successfully 
applied in Ref. |ll| to gauge theories. This cou- 
pling is defined by truncating the weak order ex- 
pansion of the plaquette (measured in the Monte 
Carlo simulation) after the first order term, and 
inverting the resulting relation. As a check, we 
investigate an alternative j3 E scheme || , defined 
by truncating the expansion after the second or- 
der term, instead. The improvement can be seen 
from Table where A^ corresponds to the /3e, 

(2) (2) 

and A^ to the (3 E scheme. 

Extrapolations to the continuum limit were 



done by fitting the data to the leading order ex- 
pectation A£ x (a) = A z ; 1 (0) + C/(Vo : ln(Z?aVo : )). 
There is a considerable bandwidth in the extrap- 
olations, depending on the number of points in- 
cluded in the fits. The effective schemes help to 
decrease this uncertainty and teach us that linear 
extrapolations are misleading. Our best guesti- 
mate for A L (0) is y/a = 50.8J^Ai. 

An interesting program would be to take other 
nonperturbative quantities (like small Wilson 
loops) that can be measured accurately on small 
lattices, and use them to define different effective 
couplings. Since one can translate, at a given (3, 
one scheme into another, self consistent extrapo- 
lations (and interpolations) become possible. By 
exploiting all available nonperturbative informa- 
tion, one can in this way reconstruct the under- 
lying /3-function beyond the two-loop approxima- 
tion. 

Running coupling: So far we have computed 
a mass (y/cr) as function of a coupling. In the 
following, we reverse this procedure and deter- 
mine the coupling a qq -(R). This is achieved [lj] 
by numerically differentiating the (lattice cor- 
rected) potential, and multiplying the resulting 
interquark force with jR 2 . The results (Figure 3) 
are in fairly good agreement with ajrg(a^ 1 ) 
( a Ms(«) = a id r ) + 0.05a^(r-) + • • • with q = 
1/r), calculated from the scaling of the string ten- 
sion in the manner proposed in Ref. JL0|. 

We fit our data to the expectation a~ q (R) = 
4tt (b Q \n(RaA R )- 2 + fei/fcolnln^aA^)- 2 ). Ex- 
ploiting the relation An = 30.19Al |l3| we ar- 
rive at y/a — 53.7(2. 1)A£ which is in agreement 
with the value extrapolated from the scaling of 
the string tension. We find the data to be in 
good agreement with the perturbative formula for 
q > 5y/a rs 2 GeV. At small energies (large R) 
our data tends to lie above the two-loop expres- 
sion since a q q oc l/q 2 (q — > 0). 

Averaging the two independently calculated A 
parameters gives y/a = 51.6J"J;gAi or 

A Ws = 0M8±?y/<r = 246tgMeV. (2) 

The approach towards this asymptotic limit for 
different definitions of effective couplings is illus- 
trated in Figure 4. We observe that the values 
for A~[ , extracted from the running of the cou- 
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Figure 3. The running coupling a q q(R) 
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Figure 4. Violations of asymptotic scaling. 



pling p],plf| , tend to increase with the momentum 
cutoff. 

4. Conclusion 

For our results, it has been important to study 
both infrared, and ultraviolett aspects in order to 
verify the reliability of the continuum extrapola- 
tion. We might say that we have been lucky to 
get hold of asymptotia within our means. This 
is due to the discovery that the running coupling 
"constant" is well described within this theory by 



the two-loop formula down to a scale of 1-2 GeV. 
If nature continues to be nice to us it is possi- 
ble to predict experimental numbers like cts(Mz) 
or Ai-L as explained in Ref. ||. In order to re- 
move additional systematic uncertainties caused 
by the fact that experimentalists still have not 
discovered how to switch off some quark flavours 
and deliver us a quenched value for the string ten- 
sion, it is preferable to repeat this study in full 
QCD on the level of TERAFLOPS power. Mean- 
while, further improvements of lattice techniques 
are of great interest. 
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